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The Influences ot .he tolling speed, inlet pressure and inlel tangent.al velocity have 
been already investigated. But the influence of the number of flns have not been experi- 
mentally investigated yet. The purpose of this paper is to investigate the influences of the 
number of flns under the rotating speed, whirling speed, inlel pressure, inlet tangentia 
velocity on a stability of the rotor system. 

The fluid force is measured unde, the condition ofboth rotating and whirling motions 
of the rotor. The rolordynamic coelhcients, which represent the flow induced force of 
labyrinth seal by the function of velocity and displacement, are obtained. 

Test Facility 

The test facility to measure the fluid force is shown in figure 1. In this figure th 
rotor is driven by the three-phase motor® . The rotor is supported at the left and right 
sides by the eccentric bearings. The dimensions and directions of both eccentricities 
are set to be equal. Therefore, the rotor is whirled by the rotation of the eccen nc 
bearing. The eccentric bearings are driven by the three-phase motor® throng e 
timing-belt and pulley® . Since the rotor and the eccentric bearing are driven by 
different motors, the rotor can be whirled in the same direct, on(forward) and in the 
opposite direction(backward), and these rotating and whirling speeds are continuously 
changed from 200(rpm) to 2000(rpm). 

The working fluid is compressed by the comprrsso.® and transported from the Str- 
unk® into the test casing through the regulator®. The test seals are changeable by 
changing seal ® and® . The detail seal dimension m this test is shown in table . 

The rotating speed of the rotor is measured by the pulse notched on the coupling. 
The whirling speed of the rotor is measured by the pulse of the balance weig t®. 

The fluid forces generated in the labyrinth seal are measured by the pressure trans- 
ducers as shown in figure 2 and six loadcells which directly support the seal as shown m 
figure 3. By measuring the force with six loadcells it is possible to measure the lateral 
force and moment. The measured signal current is sent to the AD converter roug 
the D.C. amplifier and is analyzed. Figure 7 shows the data flow chart. 

Analysis of Measured Data 

Figure 4 shows the definition of displacement, velocity and flow induced forces. 

(a) Data analysis of the prenmre transducers 

Since the data measured by the pressure transducers are the pressure distributions 
around the rotor, we can obtain the fluid forces by integrating the pressure distributions 

around and along the rotor. 

I, is assumed that the circumferential pressure distribution is elliptic and the center 
of it is O' as shown in figure 5. Thus, the circumferential pressure distribution ,n the . 
th chamber P,(0,t) is given as; 

P,(0,t) = 7 y i + AP •cos(Oi-0 + <fc) ^ 
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where T, is the magnitude of steady pressure in the »th chamber, A P, is the magnitude 
of the pressure fluctuation in the »th chamber and is the phase of the fluid force relative 
to the rotor displacement. The x and y direction components of the fluid force in the tth 
chamber are obtained by using Eq.(l) and considering the condition that the eccentric 

direction of the rotor agrees with x direction, that is, t = 0 , the following equations 
are obtained. 

2 r 

P t {9, t) ■ Ptg.Rs- cos Ode 

— — t • Rs ■ Ptg • A P, cos <j>, ^2) 

, 2 * 

Py* — ~ \ P%{9, t) • Ptg ■ Rs • sin 6d0 

Jo 

= — tt ■ Rs • Ptg ■ A Pi sin <f> t ^ 

where Rs is the rotor radius, Ptg is the seal pitch. Then, since the overall fluid forces 
generated in the labyrinth seal are calculated by adding the fluid force in each chamber, 
the following equations are obtained. 



Fx = -w ■ Rs ■ Ptg -J^AP, cos & 

»=1 

(4) 

n 

Py = — t • Rs • Ptg • ^ A P, sin <f>j 

i'=i 

(5) 


where n is the number of fins. 

(b) Data analysis of the loadcells 

The fluid forces are calculated by dividing the each forces f } measured by six loadcells 
into x,y components and adding the divided force components in the each direction. 
That is to say, the fluid forces F x, Fy is calculated by the following equations; 


Fx = f E />*«»(«>) ( 6 ) 

Fy = fE/i •““(“>) (7) 

where a 3 is each loadcelPs setting angle. 

Then, using these forces measured by loadcells, moment can be calculated. 

Figure 6 shows the definition of the moment generated in the labyrinth seal The 
moment are defined by Kanemori[9] as following equations; 

Mx = l Pyin-F yout ).h. 

it 

= ~{Fx ln -F Xout ).h 

it 


My 
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( 8 ) 

(9) 



where M x y My is the moment with regard to the x and y axes, respectively, Fx xn} Fy tn 
are forces in x and y directions measured by inlet side loadcells, Fx^t, Fy^t are forces 
in x and y directions measured by outlet side loadcells and L & is the distance between 
the inlet side and the outlet side loadcells. 

The rotordynamic coefficients are defined by the following equations. 

-{£}- 

where ( x , y ) are the rotor displacement, (F x ,F y ) are the x and y direction components 
of the reaction force acting on the rotor, and K xy} IC yxi K yy } and (Cj^, Cp, 
are the stiffness and damping, respectively. 

For small motion about a centered position, Eq.(10) is rewritten simpler as follows; 




K xx 

~K*v 



For centered position, the rotor displacement ( x,y ) are described as; 


x = ecos(ftt) 
y — esin(flt) 


( 11 ) 


( 12 ) 

(13) 


where e is the eccentricity of the rotor. 

Substituting Eq.(12), Eq.(13) into Eq.(ll) and considering the condition that the 
eccentric direction of the rotor agrees with x direction, that is, t = 0 ,the folllowing 
equations are obtained. 


- — - C xy Q + K xx 
e 

Fy 


= C X M - K 


xy 


F = \J Fx - 2 + Fy 2 , 4> = tan 1 


(14) 

( 15 ) 

( 16 ) 


where F is the magnitude of the fluid force. 

As we obtain the fluid forces for 10 whirling speed U in experiments, we can determine 
the rotordynamic coefficients by the least squares method. 


Experimental Results 

The Tange of parameter in this test is shown in table 2. The ratio of inlet tangential 
velocity to rotor surface velocity (inlet tangential velocity ratio) is shown table 3. Figure 
8 shows the comparison of the fluid force versus whirling/ rotating ratio between the 
loadcell data and the pressure transducer data. This figure also shows the magnitude 
and the phase of the fluid force measured by two ways almost agree. After this, the data 
put in this paper are the one measured by the loadcells. 
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number of fins is constant. 
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Rotordynamic Coefficients Results 
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Figure 18 show that an effect of the number of fins on the rotordynamic coefficients 
changing the inlet pressure. 

The figure 18(a) shows that the magnitude of direct stiffness increases with increasing 
inlet pressure This figure also shows that direct stiffness changes positive to negative 
with increasing the number of fins. 

The figure 18(b) shows that an increase of inlet pressure makes the magnitude of 
cross-coupled stiffness large. This figure also shows that the magnitude of cross-coupled 
stiffness decreases with increasing the number of fin. 

The figure 18(c) shows that direct damping increases with increasing inlet pressure. 
The figure 18(d) shows that an increase of inlet pressure has no influence on cross- 
coupled damping. The magnitude of cross-coupled damping increases with increasing 
the number of fins. 

Figure 19 show that an effect of the number of fins on the rotordynamic coefficients 
changing the inlet tangential velocity. 

The figure 19(a) shows that a change of inlet tangential velocity has no influence 
on direct stifness. The figures also show that direct stiffness changes from positive to 

negative with increasing the number of fins. 

The figure 19(b) shows that cross-coupled stiffness changes from negative to posi- 
tive with increasing inlet tangential velocity, and the quantity of cross-coupled stiffness 
fluctuation in this time is greater with decreasing the number of fins. 

The figure 19(c) shows that direct damping is constant with increasing inlet tangential 

velocity. 

The figure 19(d) shows that cross-coupled damping changes from negative to posi- 
tive with increasing inlet tangential velocity and the quantity of cross-coupled damping 
fluctuation in this time is greater with decreasing the number of fins. 

Figure 20 shows that an effect of the number of fins on ’whirl frequency ratio* changing 
the inlet tangential velocity. 

Whirl frequency ratio / is defined by Childs[4] as following equation; 

A'jy 

; C XX Q 

where fi is rotor whirling frequency, K xy /C xx Q is the ratio of the destabilizing influence 
of the cross-coupled stiffness and the stabilizing influence of direct damping. From 
viewpoint of stability, a less value of whirl frequency ratio is desireble. 

The figure shows that 11 fin’s labyrinth seal has a greater whirl ratio than other seals 
in giving inlet tangential velocity. And this agTees with Scharrer’s theoretical results[6]. 

Conclusion 

In this paper, the fluid force, moment and the rotordynamic coefficients are experi- 
mentally obtained. The conclusions are as follows; 

1. The labyrinth seal which has less seal fins is sensitive to increase of inlet tangential 
velocity. 
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2. The labyrinth seal which has less seal fins is not always stable. 

3. Tangential fluid force increases with increasing inlet pressure and inlet tangential 
velocity. 

4. Radial fluid force increases the magnitude with increasing inlet pressure. 

5. Moment with regard to the * axis is sensitive to the change of inlet tangential 
velocity. 

6. Moment with regard to the y axis increases the magnitude with increasing the 
number of fins. 


7. Cross-coupled damping and stiffness increases from negative to positive with in- 
creasing inlet tangential velocity. 

8. Direct damping is constant with increasing inlet tangential velocity. 
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Table. 1 Seal Dimension in this study 


Seal Clearance 

mm 

0.5 

Fin Height 

mm 

4.0 

Seal Pitch 

mm 

6.0 

Fin Thickness 

mm 

1.0 



80.0 

Number of Fin 

5,11,21 

Fin Side 

On Stator 


Table.2 Experiment Condition 


Rotating Speed w(rpm) 

1000 , 1500 , 2000 

Whirling Speed fi(rpm) 

± ( 0.2 , 0.4 , 0.6 , 0.8 , 1.0 ) xu 

Inlet Pressure P,„(MPa) 

0.2 , 0.3 , 0.4 

Inlet Tangential Velocity w tn (m/s) 

-34.0 , 0 , 34.0 


Table 3 Preswirl Velocity Ratio 



-34.0(m/s) 

0(m/s) 

34.0(m/s) 

lOOO(rpm) 

-4.06 

0.0 

4.06 

1500(rpm) 

-2.70 

0.0 

2.70 

2000(rpm) 

-2.03 

0.0 

2.03 
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Fig. I Test facility 
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Pressure transducer 



(a) Detail ol hole for pressure transducer 


(b) Cross section of 

test facility setting pressure transducers 



Fig. 2 Setting of pressure transducer 



Air Exit 



Detail of loadcells supporting labyrinth seal 


Fig. 3 Setting of loadcells 
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Fig. 4 Definition of displacement, Fig. 5 

velocity and flow induced force 


Pressure distribution around the rotor 


Y Inlet Side 



Fig. 6 Definition of flow induced force and moment 
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Pressure 

Transducers 


Loadcells 


Gap Sensor 
(for whirl) 


Gap Sensor 
(for eolation) 



Fig. 7 Data, flow chart 
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Comparison of loadcells data and pressure transducers data 

uj - 2000 (rpm) , P tn = 0.1 (MPa) , w in = 0 (m/s) 

L : Loadcells data 

P : Pressure transducers data 
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Whirl 
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(c) 5 fins W W,n = 34 0 ( m / s ) 

Fluid force versus whirling/rotating ratio Fig. 12 Fluid force versus wlnrling/rotating ratio 
(Parameter: Inlet tangential velocity) (Parameter: The number of fins) 

u = 2000 (rpm) , P in = 0.2 (MPa) " = 1000 (rpm) , P m = 0.3 (MPa) 












Circumferential Velocity [m/e] 



(a) 21 fins 






Fig. 13 Circumferential velocity versus seal length 
(Parameter: fin side) 

" = 2000 (rpm) , P in = 0.3 (MPa) 


■ Teeth on stator 

— Teeth on rotor 

io, : Rotor surface speed 
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(b) 

11 fins 

(b) 11 fins 




HMr I ing/Ratatlng 

(c) 5 fins 

Moment versus whirling/rotating ratio 
(Parameter: Rotating speed) 

= 0.3 (MPa) , Wi n = 0 (m/s) 

1 : lOOO(rpm) 

2 : 1500(rpm) 

3 : 2000(rpm) 


HH 1 1 1 t ng/Re tiling 

(c) 5 fins 

Fig. 15 Moment versus whirling/rotating ratio 
(Parameter: Inlet pressure) 
w = 1500 (rpm) , w in = 0 (m/s) 

1 : 0.2 (MPa) 

2 : 0.3 (MPa) 
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Whirl ing/Rot at ing 


(b) 11 fins 


Whirl Ing/Rat at tng 

(b) Win = 0.0 


2.5 h 



-5 I 

-1 -.5 • 

Whirl Ing/Rot at t ng 

(c) 5 fins 

Fig. 16 Moment versus wliirling/rotating ratio 
(Parameter: Inlet tangential velocity) 
w =1000 (rpm) , P xn = 0.2 (MPa) 


Whir I Ing/Rat at Ing 

(c) w in = 34.0 (m/s) 

Fig. 17 Moment versus whirling/rotating ratio 
(Parameter: The number of fins) 
w = 1500 (rpm) , P in = 0.3 (MPa) 


1 : -34.0 (m/s) 

2 
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1 : 21 fins 

n . r 







Cxx N's/m -SL* . . Kxx kN/m 


200 
100 
0 

100 
200 

Direct stiffness versus inlet pressure (b) Cross-coupled stiffness versus inlet pressure 

300 
200 
100 
0 

-100 
-200 
-300 

(c) Direct damping versus inlet pressure (d) Cross-coupled damping versus inlet pressure 

Fig. 18 Rotordynamic coefficients versus inlet pressure 
(Parameter: The number of fins) 
w = 2000 (rpm) , w in = -34.0 (m/s) 

1 : 21 fins 

2 : 11 fins 

3 : 5 fins 
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(a) Direct stiffness versus inlet tangential velocity 



(b) Cross-coupled stiffness versus inlet tangential velocity 




(c) Direct damping versus inlet tangential velocity (d) Cross-coupled damping versus 

inlet tangential velocity 

Fig. 19 Rotordynamic coefficients versus inlet tangential velocity 
(Parameter: The number of fins) j , 21 fins 

w = 2000 (rpm) , P in = 0.2 (MPa) 2 . n fins 

3 : 5 fins 
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Swirl Speed m/'s 
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: 21 fins 
: 11 fins 

: 5 fins 


500 (rpm) 
1000 (rpm) 


Fig. 20 Whirl frequency ratio versus inlet tangential velocity 
(Parameter: The number of fins) 
w = 2000 (rpm) , Pi n — 0.2 (MPa) 
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